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ABSTRACT 

Recent theoretical and observational studies have shown that ashes from thermonuclear burning may 
be ejected during radius-expansion bursts, giving rise to photoionisation edges in the X-ray spectra. 
We report a search for such features in Chandra spectra observed from the low-mass X-ray binary 
4U 1728-34. We analysed the spectra from four radius-expansion bursts detected in 2006 July, and two 
in 2002 March, but found no evidence for discrete features. We estimate upper limits for the equivalent 
widths of edges of a few hundred eV, which for the moderate temperatures observed during the bursts, 
are comparable with the predictions. During the 2006 July observation 4U 1728—34 exhibited weak, 
unusually frequent bursts (separated by < 2 hr in some cases), with profiles and a- values characteristic 
of hydrogen-poor fuel. Recurrence times as short as those measured are insufficient to exhaust the 
accreted hydrogen at solar composition, suggesting that the source accretes hydrogen deficient fuel, 
for example from an evolved donor. The detection for the first time of a 10.77 min periodic signal 
in the persistent intensity, perhaps arising from orbital modulation, supports this explanation, and 
suggests that this system is an ultracompact binary similar to 4U 1820—30. 

Subject headings: stars: neutron — X-rays: binaries — X-rays: bursts — X-rays: individual(4U 
1728-34) 



1. INTRODUCTION 

One of the highest priorities for observational stud- 
ies of neutron stars is measurement of the mass and 
radius, sufficient to constrain the uncertain equation of 
state (EOS; e.g. lLattimer fe Prakashll2007l) . A promis- 
ing avenue to achieve such measurements is to detect 
surface spectral features from accreting neutron stars, 
which may preferentially show such features during ther- 
monuclear (type-I) bursts, caused by unstable ignition of 
accreted H/He on the surface of the neutron star. Iden- 
tifying such features allows measurement of the surface 
rcdshift, and hence the compactness (M/R ratio), which 
in turn allows constraints to be placed on the EOS. In 
recent years, only one cla im for such features has been 
made (jCottam et al.ll2002T ). and subseque nt efforts have 
not b een able to confirm this result (e.g. iCottam et al.1 
20080- The majority of observational efforts to date 
have focussed on sources which show frequent bursts, 
which tend to be faint. Another class of bursts are sig- 
nificantly brighter, since the flux during the burst rise 
exceeds the Eddington limit (e.g iLewin at 
which point the radiation pressure of the burst luminos- 
ity exceeds the local gravity. Recent theoretical work 
suggests that such bursts could drive a wind contain- 
ing heavy-element ashes from the burning, which would 
imprint absorption fea tures on the X-ray s pectrum at 
the peak of the burst (|Weinberg et al.l 120061 ). Observa- 
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tional support for this hypothesis has been found from 
the low-resolution spectra of extremely intense "super- 
expansion" bursts. These bursts are the most intense of 
those thought to be powered by H/He, in contrast to the 
even longer, infrequent "superbursts" , w hich are thought 
to be powered instead by carbon (e.g. iCumming et al.l 
2006). In superexpansion bursts observed by the Rossi 
X-ray Timing Explorer (RXTE), lin 't Zand fc Weinberg! 
(2010) found highly significant residuals for fits with a 
blackbody model. The residuals could be explained by 
absorption edge features, with energies consistent with 
highly-ionised iron-peak elements and depths indicative 
of abundances > 100 times solar. However, the rela- 
tively poor spectral resolution prevents identification of 
the features. 

Here we describe an attempt to detect and measure the 
X-ray spectrum at the peak of radius-expansion bursts 
using the Chandra High energy Transmission Grating 
Spectrometer (HETGS). Few sources reliably show fre- 
quent, radius expansion bursts, and the key challenges for 
this search are the identification of a target and the ob- 
servation triggering. One candidate is the globular clus- 
ter source 4U 1820—3 which orbits in an ultracompact 
binary (P orb = 685 s iStella et all 1 19871 ) and thus likely 
accretes almost pure helium. This system only exhibits 
bursts w hen in the minimum of its ~ 180-d intensity cy- 
cle (e.g. iPriedhorskv fc Terrell U984) . and on one occa- 
sion was observed to produce 7 rad ius-expansion burst s 
separated by 193 min on average (jHaberl et alJfl987h . 
However, frequent bursts are not consistently observed, 
and recent pointed observations targeting the long-term 
minima have had limited success, for example detecting 
only four bursts i n 1 Ms of Rossi X-ray Timing Explorer 
observations (e.g. iGallowav et aLll2008l ). A better target 
is 4U 1728—34, a well-known burst source for which the 
accretion rate does not reach such extreme values, and 
is typically ~ O.lAfEdd- This rate is substantially lower 
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(by almost an order of magnitude) than the maximum 
reached by 4U 1820-30 and 4U 1728-34 commonly pro- 
duces radius-expansion bursts every « 4-5 hr. The burst 
properties are ot herwise similar to those of 4U 1820—30 
(jCumm ing 2003j), and as we show here, it is possible that 
4U 1728—34 is also an ultracompact system. 

2. OBSERVATIONS & ANALYSIS 

We observed 4 U 1728-34 between 20 06 July 17-23 
with the HETGS (jCanizares et all 120051 ) aboard Chan- 
dra. The HETGS consists of two separate grating ar- 
rays (the Medium- and High-Energy Gratings, MEG 
and HEG) which intercept about 40% of the incoming 
photons and disperse them linearly along the Advanced 
CCD Imaging Spectrometer (ACIS). These instruments 
together provide spectra in the range 1.2-31 A (10-0.4 
keV) with spectral resolution (E/AE) of up to 1000, and 
wavelength resolution of up to 0.012 A. We adopted the 
continuous clocking (CC) observing mode, for which the 
ACIS-S CCDs are read out continuously. This sacrifices 
one axis of spatial imaging but allows much higher in- 
cident photon countrates before photon pileup (multiple 
photons arriving at the same pixel in the same frame 
time) becomes a problem. Such mitigation strategies are 
necessary to deal with the highly variable incident count 
rates for a bursting source, but complicates the data re- 
duction and spectral analysis. Specifically, CC-mode ob- 
servations prevent the usual order sorting for dispersed 
spectra, which rely (in part) on both spatial axes, al- 
though the photon energy measured by the CCD allows 
unambiguous order sorting, at least for low dispersion 
orders. CC-mode observations also prevent extraction of 
the usual background spectrum. 

We reduced the Chandra data with c iao version 
4.0 (December 2007: iFruscione et all [2006] ) . using the 
calibration database (CALDB) version 3.4.2. We 
undertook spectral fitt ing with ISIS version 1.4.9-4 
(|Houck fc Denicolall2000l ). We calculated response ma- 
trices (using MKGRMf) and ancillary response files (us- 
ing FULLGARf) for each observation. We also used the 
response and ancillary files calculated for tie entire ob- 
servation to analyse tie spectra from any bursts detected 
in tlose observations. 

We also analysed arclival Chandra timed exposure 
(TE) mode observations of 4U 1728-34 made on 2002 
Marcl 4, as well as arclival XMM-Newton observations 
from 2002 October 29. We analysed tie XMM-Newton 
data using tie Science Analysis Software (SAS) ver- 
sion 7.0.0 (June 2006). We used data from tie Euro- 
pean Ploton Imaging Camera (EPIC) pn CCD cam- 
era (|Kuster et all 12002). wlicl was operated in timing 
mode. Tie EPIC pn camera is sensitive to plotons in 
tie range 0.15-15 keV, and witl an effective area peak- 
ing at around 1400 cm 2 at 2 keV. As witl CC-mode for 
Chandra, in pn timing mode tie central CCD alone is 
read out continuously, providing an effective time resolu- 
tion of 0.03 ms. Tie observations used in tlis paper are 
summarised in Table [T] 

3. RESULTS 

We obtained in 2006 July tlree separate pointings of 
approximately 50, 150 and 40 ks eacl, spanning almost 
6 days, and detecting four, eiglteen, and tlree bursts in 



eacl segment (respectively). We first analysed tie per- 
sistent source spectrum between tie bursts, to establish 
tie persistent X-ray intensity. Tie raw ligltcurve indi- 
cated significant variations in tie persistent count rate 
from observation to observation, and also between eacl 
inter-burst interval. Tlus, we divided eacl observation 
into segments delineated by eacl pair of bursts, exclud- 
ing tie intervals up to 150 s after eacl burst and 50 s 
before tie next burst. For tie first bursts of eacl point- 
ing, we integrated from tie beginning of tie pointing to 
50 s before tie burst. We tlen fit tie summed, first- 
order MEG and HEG spectra simultaneously in tie en- 
ergy range 2.1-8.3 A(6-1.5 keV) witl an absorbed comp- 
tonisation model, motivated by earlier broad-band spec- 
troscopic measurements (see iMisanovic et al.ll2010l lere- 
after M10, for full details of the persistent spectral fitting, 
and references theirein for additional analyses of the per- 
sistent spectrum seen by Chandra). The neutral column 
density was fixed at 2.29 x 10 22 cm -2 . 

Like many low-mass X-ray binaries (LMXBs), 
4U 1728—34 shows substantial variation in its X-ray 
spectrum and inten sity on a range of timescales, from 
days to decades (e.g. iLewin et al.lll993l ) , and the spectral 
"states" influence the burst properties. The integrated 
persistent flux during the 2006 Chandra observation was 
in the range 1.16-2.30 x 10~ 9 ergs cm -2 s _1 (1.5-6 keV), 
equivalent to at least 3-5 x 10 -9 ergs cm -2 s _1 in the 2.5- 
25 keV band, based on the extrapolated spectral model. 
These values are atypically high for the source; extensive 
earlier observations with RXTE f ound the source in thi s 
flux range only 10% of the time (|Gallowav et al.l [20081 ). 
In the absence of spectral information above 10 keV, 
it is difficult to identify the source spectral state; how- 
ever, the nearest RXTE observation, on 2006 July 2^3, 
found a characteristically "soft" spectrum (e.g. with 
electron temperature for a compTT spectral component 
of 4.6 keV), and X-ray colors indicating the source was 
in the lower " island" region of the X-ray color-color dia- 
gram (see e.g. [GaTlowav et al.112001 . Tie source flux was 
liglest during tie longest (second) pointing, consistently 
above 2 . x 1 ~ 9 ergs cm - 2 s -1 (1.5-6 keV). In contrast to 
previous observations (e.g. lD'Ai et al.l2006l ). no iron lines 
were detected in tie persistent spectrum. For a distance 
to tie source of 5.2 kpc (from tie peak flux of radius- 
expan sion bursts observed by RXTE] I Galloway et all 
120081 ) , and a bolometric correction wlicl assumes a char- 
acteristic high-state "soft" spe ctrum (with electron tem- 
perature of 3 keV, c.f. with iFalanga et all 120061 ) . the 
corresponding range of accretion rate is 4-7/d of the 
Eddington rate (« 1.8 x 10 18 gs _1 , assuming pure He 
accretion). 

3.1. The X-ray bursts 

We detected 25 bursts during the 2006 July observa- 
tions (Fig. [lj. The four bursts in the first observa- 
tion (#6568) were separated by 3.7 ± 0.3 hr on aver- 
age, while in the second and third observations (#6567, 
7371) the mean recurrence time was significantly shorter, 

6 Observation ID 92023-03-72-00 

7 Note that the fraction will be correspondingly higher for 
nonzero H-fraction, for example 6-11% of the Eddington rate 
A^Edd H = 10 18 gs — 1 , assuming H-fraction in the accreted fuel 
of X = 0.7. 
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TABLE 1 

X-RAY OBSERVATIONS OF 4U 1728- 
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Data 


Exposure 


No. of 




No. 


Date 


Instrument 


obsid 


mode 


(ks) 


bursts 


Ref. 


1 


2002 Mar 4 


Chandra 


2748 


TE 


29.4 


2 


[1] 


2 


2002 Oct 29 


XMM-Newton 


0149810101 


Timing (PN) 


28.1 





3 


2006 Jul 17-18 


Chandra 


6568 


CC 


49.5 


4 


[2, 3] 


4 


2006 Jul 18-20 


Chandra 


6567 


CC 


151.8 


18 


[2, 3] 


5 


2006 Jul 22-23 


Chandra 


7371 


CC 


39.7 


3 


[2, 3] 



References. — 1. ID'Af et al.l pOOBT ): 2. this paper; 3. IMisanovic et al.l ^UWi ) 

3.2. Search for spectral features during the bursts 



at 2.4 ± 0.5 hr. The bursts themselves were short, with 
rise times of « 1 s and lasting only rs 15 s (Fig. [2]). 
The four bursts in the first segment (observation #6568) 
were the brightest, with the bursts in the second and 
third segments reaching count rates only approximately 
2/3 the peak as in the first segment. A more detailed dis- 
cussion of the burst properties and the inferred ignition 
conditions can be found in M10. 

Two of the bursts in the first observation (#6568) had 
indications of a constant intensity at maximum, lasting 
w 1 s. Since constant luminosity in the peak is an indi- 
cation of radius-expansion, we carried out time-resolved 
blackbody spectral fits of the combined first-order HEG 
and MEG data for each set of bursts in each observation, 
to test for the characteristic variations in the blackbody 
radius and temperature expected for radius-expansion 
bursts. Following the conventional approach, we sub- 
tracted the pre-burst emission as background, and esti- 
mated the bolometric flux based on the blackbody spec- 
tral parameters. The average peak flux for the bursts in 
the three observations (#6568, 6567 and 7371) was 8 ±2, 
7 ± 2 and (5.9 ± 1.9) x 10 -8 ergscm _2 s _1 , respectively. 
We found evidence for an elevated blackbody radius, co- 
incident with a decrease in the blackbody temperature, 
only for the first four bursts (Fig. [3]). While we did not 
fully resolve the rise in blackbody radius, this may largely 
have been due to the relatively low count rate early in 
the burst coupled with the steep rise time. These varia- 
tions suggest that the first four bursts exhibited radius- 
expansion. Although the confidence interval for the mean 
fluxes for the subsequent bursts were consistent with the 
maximum reached by radius-expansion burs ts observed 
by other instruments (|Gallowav et al.l [2003D . the time- 
resolved spectral variations in those bursts gave no indi- 
cations of radius-expansion. 

We measured the ratio a (e.g. IGottwald et al.l Q.986) 
of the integrated persistent flux within each interval be- 
tween the bursts, to the integrated flux of the subsequent 
burst (from time-resolved spectral fits; see M10) and ob- 
tained values which were 190 in the mean. The rapid 
timescales of the bursts, as well as the high a-values, 
are both consistent with ignition of primarily He fuel. 
Although both the persistent flux and the burst recur- 
rence times varied significantly, we found no correlation 
between the two, as might be expected (under the usual 
assumption that the persistent flux tracks the accretion 
rate). As discussed by M10, this lack of correlation may 
result from incomplete burning, or alternatively a vary- 
ing area over which accretion takes place on the neutron 
star. 



We made a search for discrete features in the spectra 
from the bursts observed by Chandra in 2006 July. We 
first focussed on the four bursts observed in the first ob- 
servation, #6568, since only these exhibited any evidence 
for photospheric radius-expansion. We extracted spectra 
from various time intervals of each burst, and combined 
the ±1 orders to obtain a single first-order spectrum for 
each grating arm (HEG and MEG). We searched spectra 
extracted from 0-2 s (relative to the burst start time), 
and 0-4 s. The combined total counts in HEG/MEG 
first-order for each of these intervals was 926/912 and 
1740/1591. We grouped the bins by various factors to 
obtain bin widths of 0.01, 0.03, 0.05, 0.11 and 0.22 A, 
in order to search for both broad and narrow features. 
Since the counts per bin for these spectra was < 10 at 
short and long wavelengths (away from the peak sensitiv- 
ity) we adopted Gehrels weighting to determine errors on 
each bin, which better approximates Poisson uncertain- 
ties for low counts. This is preferable over the alternative 
approach, i.e. to bin adaptively to ensure some minimum 
number of counts per bin, since absorption features that 
saturate the spectrum and reduce the counts to zero will 
be missed. We fit each pair (HEG and MEG) of first- 
order spectra using an absorbed blackbody with column 
density uh fixed at 2.3 x 10 22 cm -2 (derived from fits 
covering the first 2 seconds of each burst; see also M10). 
We set a detection threshold of 3-er, taking into account 
the total number of bins (HEG+MEG) in each time in- 
terval and binning. 

We found no significant deviations from the blackbody 
(continuum) for any choice of binning, in either time in- 
terval. For the spectra extracted from the first two sec- 
onds of the bursts, which likely covers the entire radius 
expansion interval, we derive 3ct upper limits for emis- 
sion/absorption lines of equivalent width in the range 
100-180 eV. We show the spectrum extracted over this 
interval in Figure HI top panel. 

We also estimated the upper limits on detection of 
absorption edges in the spectra. We added edge com- 
ponents to the model with the edge energy fixed at 
w 3.1 A(4 keV), roughly where the observed count rate 
per bin of both HEG and MEG reached a maximum, 
in order to measure the best achievable limit with the 
available data. We found a 3er upper limit of 0.26 
(0.20) on the optical depth r for the 0-2 s (0-4 s) spec- 
tra respectively, corresponding to equivalent widths of 
280 (250) eV. We also searched spectra from after the ra- 
dius expansion episode, in case the photosphere was too 
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Fig. 1. — The 2006 July Chandra observation of 4U 1728—34, showing the 25 thermonuclear bursts detected, which were separated by 
recurrence times in the range 1.82-3.92 hr. The first-order photons dispersed by the HETGS were combined and binned at 0.5 s time 
resolution; the lightcurve for the zeroth-order (undispersed) photons was dominated by pilcup, which resulted in truncation of the peaks 
of the bursts, so these photons were excluded. 
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Fig. 2. — Burst profiles from the three intervals comprising the 
2006 July Chandra observation of 4U 1728-34. The first four 
bursts, from observation 6568, are plotted individually; wc subtract 
the pre-burst emission (dashed lines) and shift the profiles verti- 
cally for clarity. Dispersed events were binned to ensure 100 counts 
per bin. Wc also plot the combined lightcurves from second and 
third observations, 6567 and 7371, binned with 400 and 100 counts 
per bin, respectively. Note the markedly different peak intensity 
and shape of bursts 1-4 compared to those that followed. 
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Fig. 3. — Time-resolved spectroscopy of the summed signal for 
the first four bursts, detected in observation #6568. The histogram 
shows the (bolomctric) burst flux (left-hand j/-axis) in units of 
10 -9 ergcm _2 s _1 , with error bars indicating the la uncertain- 
tics. The grey ribbon shows the la limits of the blackbody radius 
(outer right-hand j/-axis) in km/diokpc- The diamonds show the 
la error region for the blackbody temperature (inner right-hand 
j/-axis) in keV. Note the elevated blackbody radius (to approxi- 
mately twice the value in the tail) during the second time-bin near 
the burst peak. This radius increase indicates photospheric radius- 
expansion in these bursts. 



hot during radius-expansion for neutral or incompletely- 
ionised species to survive. In spectra extracted from time 
window covering 2-16 and 4-16 s after the burst start, 
we obtained upper limits as before on edges of 250 and 
260 eV, respectively. 

Although the bursts observed in the second and third 
observations (#6567 and 7371) likely did not exhibit 
radius-expansion, we undertook a second search of spec- 
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Wavelength (A) 

Fig. 4. — HEG and MEG first-order spectra extracted from ther- 
monuclear bursts observed by Chandra. The top panel shows the 
spectra extracted from 0-2 s (relative to the burst start) from the 
first four bursts of observation #6568, in 2006 July, and binned at 
0.11 A. The summed ±1 order HEG and MEG are shown sepa- 
rately (black and red histograms), as well as the best-fitting black- 
body model (green and blue). The middle panel shows the spectra 
extracted from 0-8 s of all 25 bursts observed by Chandra in 2006 
July, binned at 0.03 A. Note the feature around 5.28 A, which is 
closest to our 3<r significance threshold, taking into account the 
number of spectral bins. The bottom panel shows the spectra ex- 
tracted from 0—5 s of the two bursts observed by Chandra in 2002 
March (observation #2748). 



tra extracted from a wider choice of time intervals (rel- 
ative to the start) of all 25 bursts. The combined total 
counts in HEG/MEG first-order for each of these inter- 
vals, i.e. 0-2, 0-4, 0-8 and 0-16 s, as well as 2-16, 4- 
16 and 8-16 s was 3540/3350, 7670/7180, 12880/12470, 
18680/19190, respectively. We found no broad or narrow 
features that exceeded our 3u confidence threshold for 
any choice of time interval or binning. The most signif- 
icant deviation was in the MEG spectra extracted from 
0-8 s of all 25 bursts, and binned at 0.03 A, in which 
a single bin around 5.28 A was below the model fit at 
a significance of almost 3er (Fig. 01 middle panel). No 
similar feture was seen in the HEG spectrum. 

Despite several attempts in recent years, observations 
of radius-expansion bursts by Chandra or XMM-Newton 
have eluded observers. In fact, we are aware of only one 
other confirmed radius-expansion burst, which was de- 
tected on 2002 March 4 22:42:45 (MJD 52337.94637) in 
observation #2748 (Table [1]) . T his burst was also de - 
tected by RXTE (burst #101 in iGallowav "eTaTl 120081) . 
and was found to exhibit photospheric radius expan- 
sion. This event was preceded by a similar burst at UT 
18:35:24 (MJD 52337.77459), 4.12 hr earlier. Under the 
assumption that this burst too reached the Eddington 
limit, we extracted spectra from the peak and the tail 
of each burst, and combined them to search for discrete 
spectral features. 

Unlike the observations in 2006 July, observation 
#2748 was performed in Timed Exposure (TE) mode, 
with the CCDs read out every 1.44 s. The extreme 
brightness of the source during the radius expansion 
bursts likely led to substantial pileup during the bursts, 
although this was not obvious from the dispersed burst 
spectra. We combined spectra from the HEG and MEG 
±1 order from the first five and next 14-19 seconds of the 
two bursts, and examined the spectra for any evidence 
of discrete features. Due to the small number of total 
counts in the two sets of spectra (600-800 counts total 
in each grating arm), we rebinned the HEG (MEG) by a 
factor of 20 (10) to achieve a bin size of approximately 
0.11 A. We obtained an acceptable fit using an absorbed 
blackbody model, and found no discrete features which 
exceeded our 3cr detection threshold (taking into account 
the number of spectral bins) . A representative spectrum, 
from the first five seconds of the two bursts summed, is 
shown in Fig. |4] (bottom panel). We estimate 3a up- 
per limits on the equivalent width of discrete features of 
±150 eV. 

3.3. Persistent intensity variations 

The long, uninterrupted X-ray lightcurve accumulated 
during the Chandra observation allowed a search for pe- 
riodic signals with sensitivity greater than for any other 
dataset, particularly at low energies. Persistent X-ray 
intensity variations in LMXBs are sometimes observed 
at the orbital period, which for 4U 1728—34 remains un- 
known. We concentrated on the longest (second) seg- 
ment of data, and created a lightcurve from first-order 
photons (both HEG and MEG) in the ks 6-25 A(2- 
0.5 keV) energy range (where orbital variations of neu- 
tral absorption would be expected to have the largest 
effect) at 60 s time resolution. We removed the 1-2 min 
segments including each burst, and filled the resulting 
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Fig. 5. — Power density spectrum for the longest contiguous seg- 
ment of the 2006 Chandra observation of 4U 1728-34. The 3- and 
4-<r detection thresholds (taking into account the number of bins 
in the frequency range of interest) are indicated, as is the bin with 
the maximum power, at a period of 10.77 min. 

gaps with normally-distributed values consistent with the 
local mean and standard deviation. The overall mean 
countrate for the lightcurve was (8.0 ± 0.4) count s _1 . 

The Fourier power spectrum of the resulting light curve 
exhibited a peak with Leahy power of 25-29 at 10.77 min, 
significant (taking into account the number of trials) at 
between 3 and 4a (i.e. at greater than 99.7% significance; 
Fig. [SJ . The peak power and significance vary because of 
the use of random values to fill the gaps in the lightcurve 
where the bursts were removed, but the detection was 
consistently > 3a over many trials. The signal was ap- 
proximately sinusoidal, with an RMS amplitude of 0.5%. 
We searched for the 10.77 min signal in the full-range 
lightcurve, as well as other energy bands. The Leahy 
power at the corresponding frequency in the FFTs from 
the full-range lightcurve was in the range 17-20, indicat- 
ing a f=s 2a detection. There was evidence that the signal 
power decreased with energy above 2 keV; nominal power 
was detected both in the 2-3 and 3-4 keV bands, with 
corresponding 3-a upper limits on the RMS amplitude of 
0.4 and 0.3%, respectively. A stronger signal with Leahy 
power rs 40 was detected in the 3-4 keV lightcurve at 
a period of 16.64 min, very close to the spacecraft yaw 
dither period of 1000 s, and so likely instrumental. We 
also searched a lightcurve comprising the zeroth-order 
photons in the < 2 keV band. The count rate in the 
undispersed low-energy photons was just 2.2 count s _1 
in the mean, substantially less than that of the first-order 
photons, likely due to pileup. We detected no significant 
power at the 10.77 min period, deriving an upper limit 
on the rms amplitude of 0.7%, which cannot rule out the 
detection in the dispersed photon lightcurve. 

Several other archival observations of 4U 1728—34 are 
available, and we also undertook timing analyses of some 
of these data in an attempt to confirm the 10.77 min pe- 



riod. The earlier Chandra observation on 2 002 March 
4 (ob servation ID #2748, Table Q] see also ID'Ai et al.l 
2006) found the source at a significantly lower intensity, 
and the first-order count rate below 2 keV was around 
5.3 count s _1 . We computed a 60-s lightcurve omitting 
the intervals surrounding the two bursts (see £ I3.2[) . and 
filled the resulting gaps with random values as with the 
2006 data. No signal was detected at the candidate pe- 
riod. With the observation only lasting 9.1 hr, the corre- 
sponding upper limit on the candidate signal amplitude 
was 1.7% (3a). An XMM-Newton observation on 2002 
October 29 (observation ID #0149810101, Tabled]) was 
similarly substantially shorter than the 2006 Chandra ob- 
servation, at only 7.5 hr. We calculated a lightcurve by 
accumulating photons with energies < 2 keV in 60-s bins, 
with a mean rate of 9.4 count s _1 . We found no signal 
at the candidate period, and derive a 3a upper limit on 
the rms amplitude of 0.8%. As with the earlier Chandra 
data, this is above the detected amplitude of the signal 
in the 2006 observations. 

3.4. The X-ray position and optical counterpart 

The detection of a short-period signal in the persistent 
X-ray emission during the 2006 observation led us to re- 
exami ne the optical counterpart proposed bv lMarti et al.l 
(1998). The infrared magnitudes combined with the 
measured column density lead those authors to suggest a 
counterpart "less luminous and less massive than a mid- 
dle B / early A V star or a F-K III star" , taking into ac- 
count the substantial contribution of the accretion disk. 
Here we investigate whether the observations are consis- 
tent with an ultracompact counterpart similar to that of 
the similar 10-min binary 4U 1820-30. 

The counterpart to that sy stem was de t ected in the 
B and F140W UV bands by King et al.l (fl99l . The 
ratio of the absorption-corrected fluxes in the two bands 
was approximately (AA/A) -3 2 , so that the spectral slope 
was close to the Rayleigh- Jeans value. We extrapo- 
lated the observed spectrum to the J and K bands, 
and corrected for the different distances of 7.6 kpc for 
4U 1820-30 and 5.7 kpc for 4U 1728-34 (assuming that 
radius-expansion b ursts reach the Eddin gton limit for a 
helium atmosphere: iGallowav et al1[20 03). The expected 
magnitudes at the distance of 4U 1728—34 are J = 17.5 
and K = 17.3, which are substantially different to the 
ab sorption-corrected values of 19.6 and 15.1 measured 
by iMartf et al.l (|1998l ) . On the other hand, if we extrap- 
olate using a spectral slope of 4, we expect J = 19.4 
and K = 19.7, which at least matches the J-band lumi- 
nosity for the proposed counterpart. We conclude that 
a counterpart similar to 4U 1820—30 may give a con- 
sistent J-band luminosity to the proposed 4U 1728—34 
counterpart, but the measured if -band brightness is too 
high. This may indicate an additional IR component as 
has been inferred to e xist in other compact binaries (e.g. 
IGreenhill et~aT1[200l . 

Alternatively, we consider the possibility that the pro- 
posed counterpart is unrelated and positionally coinci- 
dent with 4U 1728—34 by chance. First we determined 
an accurate X-ray position from the Chandra observa- 
tion #2748 (Table [J). The zeroth-order image for this 
bright source was significantly piled up, leaving a ring of 
emission rather than the point source usual for fainter 
sources. The standard source detection algorithms can- 



Radius-expansion burst spectra from 4U 1728—34 



7 



not provide a precise position for such images, so instead 
we fit the position of the dispersed HEG and MEG spec- 
tra, as well as the readout trace (an artifact arising from 
the shuffling of charge along the CCD illuminated by 
a bright point source) and computed their intersection. 
We first checked for, and corrected, an aspect offset of 
-C/46 in R.A. and -l'/03 in DecQ We determined a 
position for 4U 1728-34 of R.A. = 17 h 31 m 57!69, decl. 
= -33°50'01!3 (J2000.0), with an estimated 90% confi- 
dence uncertainty of 0'.'6 (corresponding to the system- 
atic pointing uncertainty for Chandra). We n ote that this 
posit ion is within 1'.'3 of the radio source (|Marti et al.l 
1998), confirming the positional coincidence of the radio 
and X-ray sources. 

Finally, we estimated the field density of stars with 
K < 15 based on the IR star count model of 
iNakaiima et al.l (|200Cf) , We calculate a total of 280 ob- 
jects per arcmin 2 towards 4U 1728—34, which suggests 
a reasonable likelihood of chance positional alignment 
within an arcsec. Thus, we cannot place any significant 
constraints on the nature of 4U 1728—34 based on the 
properties of the proposed counterpart. 

4. DISCUSSION & CONCLUSIONS 

The 2002 and 2006 Chandra observations represent 
the only high-spectral resolution observations of radius- 
expansion bursts from any source to date. Only six such 
bursts have been observed in total, and the data from 
two of those bursts (from observation #2748) likely suf- 
fer from pileup as well as having too poor time resolution 
to resolve the radius-expansion episode. Nevertheless, 
this work represents the first such probe into a poorly- 
studied burst spectral regime. We detected no evidence 
for discrete features in the spectra extracted from during 
the radius-expansion episodes, and derive upper limits 
for narrow features of approximately 150 eV. 

The super-Eddington flux that drives the radius expan- 
sion during these bursts transforms nearly all the excess 
luminosity into kinetic and potential energy of the ex- 
tended atmosphere. The radiation-driven winds which 
result are capable of expelling > 5% of the accrete d 
mass (e.g. lHanawa fc Sugimotolll982HNobili et al.lll994h . 
If a fraction of the heavy elements produced during burn- 
ing are expelled in a wind and exposed at the surface, we 
may expect to observe photoionisation edges from these 
elements. 

However, we found no evidence for photoionization 
edges, and derive typical upper limits of 250-300 eV. 
The strength of edge features depends upon the elemen- 
tal abundances in the burning products, the tempera- 
ture and degree of expansion of the photosphere (that 
sets the degree of ionisation) and the amount of ma- 
terial ejected (that determin es the column depth; e.g. 
lin 't Zand fc Weinberg! 120101 ). In order to estimate the 
likely equivalent widths, we adopted model predictions 
from [Weinberg et al.1 ((2006) . We assumed abundances 
in the bur ning products for the pure-He model (see e.g. 
Fig. 13 of ( W einberg et al.f2006[ )). since the burst profiles 
and other evidence presented here suggest pure He accre- 
tion. The abundance distributions may vary if there is 

8 This correction was required since wc 

used data processed prior to 2004; see 

|\protect|fcttp: / /cxc.harvard.edu /cal /ASPECT /fix.offsetl 



hydrogen present in the ignition layer. We note that the 
rp-process burning that would occur if H is present (e.g. 
ISchatz et aH ^OO!) occurs after the radius expansion, so 
heavy burning products will not be ejected in the wind, 
unless they get dredged up from the ashes of previous 
bursts, lin 't Zand fc Weinberg! (|2010( ) found evidence for 
absorption edges with optical depths of r = 0.2-3 in su- 
perexpansion bursts, in which the photospheric radius in- 
creases to > 1000 km. The expansion in the bursts from 
4U 1728—34 was much more modest, at around 20 km 
(a factor of « 2 above the NS surface; Fig. [3]), and it is 
possible that this distinction is the reason for the non- 
detection in our data. However, we note that the upper 
limits on optical depth that we derive, of r ss 0.2-0.3, 
are consistent with the measured values for the weaker 
edges in the superexpansion bursts. Thus, for moderate 
expansion bursts such as those seen from 4U 1728—34, 
we can rule out the strongest edges detected in the su- 
perexpansion bursts, but not the weakest. 

The predicted column density and equivalent widths 
for selected elements are shown in Fig. [6] For com- 
parison, we calculated the effective temperature in the 
neutron-star photosphere, adopting a corr ection to the 
obser ved (color) temperature of 2 (see e.g. iMadei et al.l 
l2004f ) and a surface gravitational redshift of 1.31 (equiv- 
alent to a 1.41 Mq neutron star with R = 10 km), that 
have the net result of reducing the observed (color) tem- 
perature by a factor of 1.5 to give the effective temper- 
ature in the neutron star frame. We calculated the ef- 
fective temperatures both at the time of maximum ex- 
pansion of the photosphere, and also at the point where 
the photosphere is assumed to have rejoined the neutron 
star surface ("touchdown"). The detectability of edge 
features is strongly dependent upon the effective tem- 
perature in the derived range, between 1 and 1.8 keV. 
Interestingly, the predicted equivalent width for S and 
Si at an effective temperature of ~ 1 keV are a few hun- 
dred to a thousand eV, comparable to the upper limits on 
detection derived in the present spectra. Doppler broad- 
ening due to rotation may reduce the detectability of the 
edges in moderate radius expansion bursts. If angular 
momentum in the expanding photosphere is conserved, 
the transverse velocity of the expanded shell in the bursts 
from 4U 1728—34 would be reduced by a factor of two, 
such that the degree of Doppler broadening would be 
equivalent to a source rotating at ~ is sp i n /2 = 181 Hz. 
The reduction in detectability for an intrinsically broad 
feature such as an edge will not be so severe as for spec- 
tral lines (e.g. lChang et al.|[2006l) . although it may be sig- 
nificant. We hope that future observations, that will con- 
tribute additional photons to the already-accumulated 
spectrum emitted at the peak of radius-expansion bursts, 
will allow us to make more sensitive searches for spectral 
features. 

The properties of the bursts observed from 4U 1728—34 
by Chandra suggest that the fuel contains little or no 
hydrogen. During the 2006 July observation the source 
entered a state of unusually high accretion rate, so that 
(from the second observation, #6567 onwards) the burst 
ignition conditions were reached before sufficient fuel to 
power a radius-expansion burst had accumulated. The 
occurrence of bursts with characteristic He-rich profiles 
and high a- values at such short recurrence times are dif- 
ficult to understand based on standard burst theory. For 
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Fig. 6. — Column density (top panel) and edge equivalent width 
(bottom panel) for representative elements predicted for model 
HcO.l (pure He accretion at 0.1 M E dd) of lWeinbere et al.l 112009) . 
as a function of effective temperature. The estimated effective 
temperature at the time of maximum radius expansion and subse- 
quent touchdown during the first four bursts observed by Chandra 
in June 2006 are indicated (vertical dashed lines). The upper limit 
on edge depth derived for the summed spectra from four radius- 
expansion bursts from 4U 1728—34 is indicated. 

a neutron star accreting mixed H/He at around 10% 
the Eddington rate (within the range in which we ob- 
served 4U 1728-34), H will burn stably between the 
bursts, reducing the fraction in the fuel at ignition. How- 
ever, the burst recurrence time observed in 2006 July 
is insufficient to exhaust the H for solar abundances. 
Steady (hot CNO) H-burning is mediated by /3-decays 
which limit the rate, so that the accreted hydrogen 
will only be exhausted at the ignition depth in a time 
t = 11('Xo/0.7)(Zcno/0.02)~ 1 hr, where X is the hy- 
drogen fraction in the accreted material and Zcno is the 
mass fraction of CNO nuclei (0.7 and 0.02 respectively, 
for solar composition). For burst recurrence times of 
w 2 hr there is insufficient time between bursts to re- 
duce the H-fraction in the fuel significantly below the 
accreted level prior to burst ignition. If the source were 
accreting hydrogen at approximately solar composition, 
the expected burst profiles would be much longer, with 
w 5 s rises and tails out to 100 s, similar to what has 
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been observed for GS 1826-24 (IGallowav et al.ll20"ol . In 
those bursts the longer timescale comes from the much 
slower progression of /3-decay mediated rp-process H- 
burning. We suggest three possible explanations for the 
burst properties. First, that the source accretes mate- 
rial which is unusually rich in CNO nuclei. With Zcno 
approximately 5 times in excess of the solar value, the 
time to exhaust all the H would be of order the burst 
recurrence time observed during the 2006 Chandra ob- 
servation. There is no observational support for this sug- 
gestion, as the optical counterpart for the source has not 
been confirmed, so the type of the mass donor is not 
known. Second, it is possible that shear-mediated mix- 
ing allows CNO nuclei from the ashes of previous bursts 
to be dredged up into the freshly accreted material, as 
has been predic t ed for slowly rotating neutron stars by 
iPiro fc Bildsten] (12007ft. 4U 1728-34 is k nown to be ro- 
tating at 363 Hz (jStrohmaver et aflll996| ) , which is fairly 
slow for accreting neutron-stars in LMXBs, but the full 
effects of such mixing in sources accreting hydrogen as 
well as helium have not yet been fully investigated. Ei- 
ther of these possibilities leave open the chance of long 
H-rich bursts occurring occasionally from the source, al- 
though these have not yet been observed. Continued ob- 
servations, as well as an archival search underway with 
data from multiple instruments, will test these hypothe- 
ses. 

Third, it remains possible that the source accretes H- 
poor material, from an evolved (e.g. white dwarf) donor. 
Such systems are referred to as ultracompact binaries, 
due to the ir typically very short (< 80 min) orbital pe- 
riods (e.g. lin 't Zand et al.ll2007h . The 10.77 min modu- 
lation in the intensity of low-energy X-rays detected for 
the first time by Chandra from 4U 1728—34 supports 
this possibility. If this modulation was orbital, the sys- 
tem would be essentially a twin of 4U 1820—30. Like 
that source, 4U 1728—34 will be detectable with the 
ESA/NASA La ser Interferometer Space Antenna (LISA; 
Nelemans 2008). Existing measurements of the candi- 
date IR counterpart cannot significantly constrain this 
possibility, and deeper infrared observations are required 
to unambiguously confirm the ultracompact nature of 
this source, via detection of a corresponding modulation 
in the counterpart intensity. 
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